Hydrodynamic correlations in the translocation of biopolymer through a nanopore: 

theory and multiscale simulations 
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We investigate the process of biopolymer translocation through a narrow pore using a multiscale 
approach which explicitly accounts for the hydrodynamic interactions of the molecule with the 
surrounding solvent. The simulations confirm that the coupling of the correlated molecular motion 
to hydrodynamics results in significant acceleration of the translocation process. Based on these 
results, we construct a phenomenological model which incorporates the statistical and dynamical 
features of the translocation process and predicts a power law dependence of the translocation time 
on the polymer length with an exponent a ~ 1.2. The actual value of the exponent from the 
simulations is a = 1.28 ± 0.01, which is in excellent agreement with experimental measurements 
of DNA translocation through a nanopore, and is not sensitive to the choice of parameters in 
the simulation. The mechanism behind the emergence of such a robust exponent is related to the 
interplay between the longitudinal and transversal dynamics of both translocated and untranslocated 
segments. The connection to the macroscopic picture involves separating the contributions from the 
blob shrinking and shifting processes, which are both essential to the translocation dynamics. 



I. INTRODUCTION 

Translocation of biopolymers, such as DNA and RNA, 
plays a vital role in many important biological processes, 
such as viral infection by phages, inter-bacterial DNA 
transduction or gene therapy [l[. This has motivated a 
number of in vitro experimental studies, aimed at ex- 
ploring the translocation process through protein chan- 
nels across cellular membranes 0; 13 > or through micro- 
fabricated channels In particular, recent experimen- 
tal work has focussed on the possibility of fast DNA- 
sequencing by "reading-off" the DNA bases while track- 
ing its motion through nanopores under the effect of a 
localized electric field [1, [6| . 

The translocation of biopolymers is a complex phe- 
nomenon involving competition between many-body 
atom-atom interactions, fluid-atom hydrodynamic cou- 
pling, as well as the interaction of the polymer with 
wall molecules in the nanopore. Although some uni- 
versal features of the translocation process can be an- 
alyzed by means of suitably simplified statistical models 
[t,, 8, 9, 10], and non- hydro dynamic coarse-grained or mi- 
croscopic models [ll], [l^, EM [13 j ^ quantitative descrip- 
tion of this complex phenomenon calls for realistic, state- 
of-the-art computational modeling. Work along these 
lines has been recently reported by several groups, begin- 
ning with the first multiscale simulations by the present 
authors H, [IB] , followed by Langevin dynamics simula- 
tions [17] and more recently by coupled molecular-fluid 
dynamics [H, [13] ■ Specifically, Forrey and Muthuku- 
mar performed Langevin dynamics simulations, and ex- 
amined single-file as well as multi-file translocation fl^. 
Izmitli et al. used a coupled lattice Boltzmann - Molec- 
ular Dynamics scheme [T^ , as outlined in [l^, [l^ , and 
reproduced our early results, while exploring a smaller 



range of chain lengths and ensemble size. These recent 
works have provided a wealth of new computational re- 
sults and detailed insight into the problem of translo- 
cation through nanopores. In addition to these studies. 
Slater and co-workers have investigated the translocation 
process numerically by also including hydrodynamic in- 
teractions [19,]. However, these authors treat the cases 
were no external field is applied, thus comparison with 
our work cannot be made at this point. 

In this work, we report a synthesis of the simula- 
tional result into a coherent mean-field analytical model 
which captures the basic physical mechanisms behind the 
translocation process. The model is based on extracting 
the scaling behavior of translocated and untranslocated 
segments, including the anisotropy between longitudinal 
and transverse components. The analytical model pre- 
dicts a power-law scaling behavior of translocation time 
with polymer length with exponent a ~ 1.2, which is very 
close to the one found in the current {a — 1.28 ± 0.01) 
and other {a = 1.28 ± 0.03 [3]) simulations and in ex- 
periments of DNA translocation {a = 1.27 ± 0.03 

The paper is organized as follows: in Section 11 we 
provide an overview of our simulational approach. In 
Section 111 we discuss the choice of simulation parame- 
ters that make our simulations relevant to DNA translo- 
cation through nanopores and the implications for the 
implied time-scales and length-scales of the system. In 
Section IV we discuss the results of the translocation sim- 
ulations, paying particular attention to the anisotropy of 
longitudinal and transverse components of the translo- 
cating polymer in both the untranslocated and translo- 
cated segments. Section V presents the analytical mean- 
field model and its comparison to the simulations. We 
conclude in Section VI with some comments on what 
may be the limitations of the mean-field picture and a 
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summary of our results. 



II. MULTISCALE MODEL 

Our multiscale method is based on the coupHng be- 
tween constrained Molecular Dynamics (MD) for the 
polymer evolution and a lattice Boltzmann (LB) treat- 
ment of the explicit solvent dynamics [U . In con- 
trast to Brownian dynamics, the LB approach handles 
the fluid-mediated solute-solute interactions through an 
explicit representation of local collisions between the sol- 
vent and solute molecules. We will focus on the fast 
translocation regime, in which the translocation time r is 
much smaller than the typical relaxation (Zimm) time of 
the polymer towards its native (minimum energy, maxi- 
mum entropy) configuration. This regime cannot be cap- 
tured by a simple, one-dimensional Brownian model [23| 
or a Fokker-Planck representation. 

Translocation is induced by the constant electric field 
E acting along the x direction localized in a region near 
the pore. The dynamics of the beads which constitute 
the molecule are governed by the equation 
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(1) 



with Ftot i the total force on bead i. F^^i is a conserva- 
tive force describing the sum of bead-bead and bead- wall 
interactions; F^^ag,! is the dissipative drag force due to 
polymer-fluid coupling given by —mh^(vi — Ui) with 7 the 
friction coefficient and Vi, Ui the bead and fluid veloci- 
ties at the position of bead i with a mass to^; Fj- i is 
a random force on bead i with zero mean; i^^^i is the re- 
action force resulting from A'o — 1 holonomic constraints 
for molecules modelled with rigid covalent bonds, with 
TVq the number of beads in the polymer; and Fdrive,i is 
the driving force representing the effect of the external 
field E^ equal to qEg(fi) with q an effective charge, which 
acts only on beads in the pore region. The region over 
which the external field acts is described by the function 
g^fi), which is 1 for fi within this region and otherwise; 
the extent of this region is chosen to be a cube of side 
3 Ax (shown by the green shaded region in Fig. [1]), where 
Ax is the lattice spacing. 

For the bead-bead interaction, other than that between 
consecutive beads, and for the bead-wall interaction we 
choose separately a truncated Lennard-Joncs potential 
(repulsive part only [28j): 



Vir) 




iir < rcut 
if r > rcut 



In this expression, r is the bead-bead or bead-wall dis- 
tance. In both cases the potentials are truncated at the 
cut-off distance r^ut = S^^^cr. The chosen parameters a 
and e are 1.8 Ax and 10~^ AmAx/At^ for the bead- 
bead interactions and 1.5 Ax and lO^"' AmAx/At^ for 




FIG. 1: (Color online) Illustration of the interactions in the 
DNA translocation model: the beads representing the DNA 
are shown as small black dots connected by straight line seg- 
ments, the wall is represented by a plane of on-lattice points 
(small white dots separated by the lattice spacing Ax) which 
repel the DNA beads within a range of interaction. The bead- 
bead interactions are indicated by the blue spheres surround- 
ing the small black ones, defined in Eq.(l2]), except for the 
distance between consecutive beads which is fixed at b. Inter- 
actions between the beads and the wall are indicated by the 
red spheres surrounding the small white ones, defined in Eq. 
([2]). Interaction between the beads and the constant external 
field are confined over a shaded (in green) region around the 
pore (see text for details). 



the bead- wall interactions. The parameters are given in 
LB units as explained in Section III. The distance be- 
tween consecutive beads along the chain representing the 
macromolecule is set to 6 = 1.2 Ax through a constraint 
imposed by the SHAKE algorithm [23] ■ The geometry 
of the pore region is shown schematically in Fig. [1] The 
fluid is represented through lattice Boltzmann particles 
that reside on a three-dimensional cubic lattice with spac- 
ing Ax. The probability distribution /p(x, t) denotes the 
number of particles at the lattice position given by x at 
time t, and evolves in space and time toward the equilib- 
rium distribution /^^ and with relaxation frequency uj, 
as: 



fp{x + CpAt,t 



At) = f.p{x,t)-uAt{U- f;'i){x,t) 



-FpAt 



GpAt 



(3) 



The LB particles can only move on the lattice with fixed 
first- and second-neighbor speeds Cp (19 in all, for a 3D 
cubic lattice), while Fp represents thermal fluctuations 



(2) and Gp describes the polymer-fluid back reaction 



Gp{x, t) 



WpP ^ 



Fr_ 



(4) 



with Wp a set of weights normalized to unity, and (3 the 
inverse fluid temperature. Finally, -D(x) denotes the lat- 
tice cell to which the i*'' bead belongs. Details on the 
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numerical implementation of the scheme have been re- 
ported in Ref. p^ . 

Compared to other numerical methods, the present 
lattice Boltzmann-Molecular Dynamics (LB-MD) scheme 
has certain computational advantages, namely, it permits 
to take into account self-consistent hydrodynamic corre- 
lations at computational cost scaling linearly with the 
polymer length. This allowed us to simulate long chains 
in 3D over large statistical ensembles at an affordable 
computational cost. We simulate polymers of various 
sizes and as large sample realizations as our computa- 
tional resources permit; specifically we considered sizes 
oi No^20 [1000], iVo=50 [1000], TVq^IOO [500], iVo=200 
[300], iVo=300 [300], 7Vo=4G0 [200], iVo=500 [150], where 
the numbers in brackets are the sample realizations. 



III. CHOICE OF THE SIMULATION 
PARAMETERS 

We discuss next the details of the choice of parameters 
in the model so that our simulations are relevant to DNA 
translocation through nano-pores, as observed in typical 
experimental setups j^]. The simulations are performed 
in a three-dimensional box which contains the polymer 
and the fluid solvent, and has a size x Nx/2 x Nx/2 
in units of the lattice spacing Ax; we used = 80 
for polymers of size iVo < 400 beads. The separating 
wall is located in the mid-section of the x direction, at 
x/Ax = Nx/2, with a square hole of side 3Aa; at the 
center through which the polymer translocates from one 
chamber to the other. At t = the polymer resides 
entirely in the right chamber at x/ Ax > Nx/2^ with its 
one end at the pore region along the streamline {x) and 
cross-flow (y, z) directions. The LB time step is At and 
the MD time step At/h. 

In order to have a plausible representation of DNA, 
we choose the separation between consecutive beads to 
be 6 = 50 nm, the persistence length of double-stranded 
DNA. Since we defined b — 1.2 Ax, this gives for the lat- 
tice spacing Aa: = 42 nm. Since the repulsive interaction 
between the beads and the wall involved the parameter 
(Tu, = 1.5 Ax, and the pore is a square of size 3Aa;, this 
produces an effective hole of size ^ 5 nm through which 
the polymer translocates. Having set the value of Ax, 
we choose the time step so that the kinematic viscosity 
is expressed as: 



J^LB- 



(Ax)2 



At 



with the viscosity of water (10^^ m^/s) and vlb the 
numerical value of the viscosity in LB units; taking vlb — 
0.1 produces a time-step of At = 160 ps. 

The above choice for the lattice spacing and time- 
step fixes the lattice speed c = {Ax/ At) ~ 250 m/s, 
w hich is r easonably close to the solvent thermal speed 
^JkT /m^ ~ 500 m/s. It is instructive to compare the 
lattice speed, c, with the typical propagation speeds of 



the main phenomena of interest in the translocation pro- 
cess. These are: 

i) the translocation speed b{dN/dt) ~ lO^'^ m/s; 

ii) the electric drift speed {qE) / {nihpjtp) ~ 0.3 m/s; 

iii) the base-pair thermal speed ^/kT/mtp ~ 70 m/s, 
where we have used the following reference values for 
the base-pair translocation rate dN/dt ^ 10^ bp/s, the 
strength of the force due to the external field qE ~ 10 
pN, the mass of a base-pair mf,p in terms of the mass of 
the water molecule , nibp ~ 30 , ^ 600 amu, the 
drag coefficient for a base-pair ji,p = 3 x 10^'^ s~^ [2§] 
and T = 300 °K for the temperature. These order- 
of-magnitude estimates indicate that the present choice 
of space and time units is such that the corresponding 
speed, u, fulfills the numerical stability CFL (Courant- 
Friedrichs-Lewy) condition, u < c. 

In LB simulations mass units are fixed by the mass 
density of the fluid species pw Setting the LB mass 
density plb = 1-0 corresponds to having a number of 
water molecules p^Ax"^ / pLBm^ — 2 y. 10^, where 
is the density of water, and each lattice site contains a 
solvent mass Am = pyjAx^ / plb- Since we are using 
a continuum-kinetic representation of fluid flow at the 
nanoscale, a necessary condition for this representation 
to hold against statistical noise is that the above ratio 
be much greater than unity. Given the fact that we have 
19 discrete distributions, fi per cell, each of them would 
represent about 10^ water molecules, a safe value towards 
satisfying this condition. 

The prime goal of the mapping procedure is to secure 
the correct values of the major dimensionless parame- 
ters governing the physics of the translocation process. 
In particular, this regards the ratio of external drive to 
thermal forces which we will call 6, defined as 



qEh 



(5) 



The value of this quantity in actual experiments is (f> ^ 
1 — 10 3. In our simulations, we took kT = 10"-* and 
qE = 0.01 acting on one bead, mapping ~ 100 base- 
pairs. These quantities are again given in LB units. At 
a base-pair level, this means qE = 10~^, corresponding 
to 4> = 1, in satisfactory order-of-magnitude agreement 
with experiments. With this value of the driving force 
and the polymers studied in the present work, a typical 
translocation event takes place in a time interval on the 
order of lO'^ — 10** LB timesteps. 

In the simulations we define the effective mass of the 
beads to be mi, = 1. A straightforward calculation shows 
that the ratio of the effective bead mass to the mass re- 
sulting from the DNA coarse graining, that is, one bead 
representing ~ 100 base pairs, is ~ 700. On the other 
hand, the parameter relevant to momentum exchange is 
the bead friction 7. We chose 7 = 0.1 in LB units, so 
that the ratio between friction in physical units and the 
experimentally determined one {■^bp = 30 ps~^ [2^) is 
10~^. Our choice was dictated by the criterion of nu- 
merical stability 7 < 1/At = 6 x 10^ s~^. These fac- 
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FIG. 2: (Color online) Scaling of the translocation time r with 
the number of beads A'^o, in the presence (circles) or absence 
(squares) of coupling between the molecule and the solvent. 
The exponents are 1.28 ± 0.01 and 1.36 ± 0.03, respectively. 



tors taken together show that the bare particle mobihty, 
fi ~ 1/ (mb7), is a factor of 15 larger than the experimen- 
tal one. This is equivalent to an underdamped motion 
of the macromolecule, which results in smooth particle 
trajectories and allows for algorithmic stability without 
affecting the long-time behavior of the polymer. A possi- 
ble alternative would be to solve the polymer dynamics in 
overdamped (Brownian) form which also circumvents the 
small timestep issue imposed by the frictional damping 
[H, [13| • Our choice to use inertial dynamics was based 
on previous experience with similar systems, without hy- 
drodynamic interactions, which showed that the inertial 
dynamics approach has a slight advantage in stability 
with larger time-step size. 



IV. TRANSLOCATION SIMULATIONS 

We report next the results of extensive simulations of 
translocation events. In each case, the translocation time 
exhibits a statistical (not exactly gaussian) distribution 
around the most probable translocation time, in close 
analogy to experimental observations 0, [1, . Accord- 
ing to these studies, a 48.5 kbp double-stranded DNA 
takes about 2 msec to translocate, while a typical Zimm 
time is about 700 msec. Fig. [2] shows the dependence 
of the most probable translocation time r on the poly- 
mer length, which obeys a superlinear scaling relation: 
T oc Nq. The observed exponent, with the molecule- 
fluid coupling in operation, a ~ 1.28 ± 0.01 is in very 
good agreement with recent experiments of DNA translo- 
cation through a nanopore driven by an external electric 



field, where a ~ 1.27 ± 0.03 (6|. 

In the absence of coupling between the molecule and 
the solvent, the translocation process is slowed down, 
as indicated by a higher exponent, a ~ 1.36 ± 0.03. 
We propose that in the presence of polymer-solvent cou- 
pling some form of bead-bead screening takes place, the 
biopolymer never translocates in the form of a linear 
chain, as this configuration is entropically suppressed. 
The physics of the larger exponent (slower motion) in the 
absence of hydrodynamics is related to enhanced corre- 
lation effects and a wider range of polymer fluctuations. 

In Fig. [3] we show an ensemble of 100 polymer config- 
urations at three different instants, referring to the ini- 
tial, mid-point and final stages of translocation. From 
this figure it is clear that, initially, the shape of the un- 
translocated (U) segment is squeezed against the wall and 
takes the form of an oblate ellipsoid. The translocated 
(T) segment appears to be more compact and prolate. 
The tendency of the U-segment to be attracted to the 
wall is in line with the well-known "mushroom" shape 
arising from a polymer anchored at one end to a repul- 
sive wall [34I . To investigate the consequences of this 
anisotropy, we have inspected the scaling of the gyration 
tensor with the number of monomers for the U and T 
segments separately, with time, according to: 



Rijt) K [Ni{t)r--,{I ^U,T; 



(6) 



where fi =||,^ denotes the longitudinal and transverse 
components of the gyration tensor, with respect to the 
direction of translocation. This scaling is shown in Fig. 
m from which we observe that for Nu > 100 the trans- 
verse component Ru.i_ follows a dynamic scaling law with 
iyu,± — 0.6, close to the Flory-exponent of a 3-D self- 
avoiding random walk, while the longitudinal one ex- 
hibits a much weaker dependence on Nu (smaller slope) . 
The T-segment follows a similar trend, but with a trans- 
verse component scaling with i^T,± — 0.5. 

The fluid/biopolymer system can be approached as 
an extended dynamical system consisting of two com- 
ponents: this system receives energy from the exterior 
through a localized electric field acting on the polymer 
and dissipates it via interaction of both the fluid and 
the polymer with the wall. Each bead is subject to the 
following forces: (a) the localized drive Fdrive,i, (b) the 
dissipative drag Fdrag.i, (c) the pore drag Fpore,i, and 
(d) the entropic forces. At equilibrimii, the latter can be 
expressed as 



ksT 1 - 2r 
b r(l — r) 



(7) 



where r{t) = NT{t)/No is the translocation coordinate Q 
and b is the separation between two beads. Here, we have 
used an explicit dependence on r simply to show that 
entropic forces are negligible most of the time, except at 
both ends (initiation and completion) of the translocation 
process. These forces are naturally measured in terms of 
the thermal force, Fth = ksT/b. In the fast translocation 
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regime considered in this work, Fdrive/Fth > 1, yielding 
bFdrive/kBT ^ 10^ for the parameters used. Pore forces 
are neghgible in our simulations, due to the small pore 
size. For the entropic forces, Fentr/Fth ~ (1 — 2r) /(r(l — 
r)), which shows that they can be neglected except at the 
early {r{t) 0) and final (r(t) 1) stages of the process. 
With entropic forces and pore dissipation negligible, the 
forces guiding the translocation are the hydrodynamic 
drag and the drive from the external electric field. 

Hydrodynamics is expected to provide a cooperative 
background, helping to minimize frictional effects. For 
quantitative insight into this, we monitor the synergy fac- 
tor, defined as the work per unit time made by the fluid 
on both parts {I—U or T) of the polymer: 



dW. 



dt 



7( > :4"it)-Mt) 



(8) 



where brackets denote averages over different realiza- 
tions of the polymer for a given length. Positive (neg- 
ative) values of S^j^\t) indicate a cooperative (competi- 
tive) solvent, respectively. The variation of S'^\t) with 
time is linear (Fig. [IJa)), while the total rate Snit) = 

(T) (U) 

(0 + '^H (0 whole chain, is constant with 

time indicating that the work per unit time associated 
with the change of the radii of the two blobs is constant. 
The probability distribution of Snit) during transloca- 
tion (Fig. [5fb)), indicates that hydrodynamics turns the 
solvent into a cooperative environment: the distribution 
lies entirely in the positive range. This cooperative effect 
is the underlying reason for the faster translocation pro- 
cess in the presence of hydrodynamic interactions. Simi- 
larly, the work done per timestep by the electric field on 
the polymer can be defined as: 



SE{t) = 



dWs 
dt 



J2F,r 



■V^{t) 



(9) 



The average of SE{t) is positive (Fig. [5]Ja)). However, 
the negative tail of the corresponding probability distri- 
bution in Fig. HKb) indicates that there is a non-negligible 
probability to find beads moving against the electric field. 
On average, Ssit) is also constant with time (except very 
near the completion of the translocation) denoting that 
the beads traverse the pore with basically the same speed 
at all times. The average hydrodynamic work per time is 
larger than the input of the external electric field, since 
the latter only acts on a very small fraction of the beads, 
about 4 resident beads within the pore region. 



V. PHENOMENOLOGICAL MODEL 

We next present a phenomenological model for the 
scaling exponent of the translocation time in the pres- 
ence of the solvent. We first consider Eq.([T|) and multi- 
ply both sides of the equation with the bead velocity Vi. 




FIG. 3: 3-D view of an ensemble of 100 polymers with A'^o — 
300 at different stages of the process: A, B and C correspond 
to the initial, mid-point, and final translocation times. 



Summing up over all beads, using the fact that the ve- 
locities are uncorrelated with the random forces F^^i and 
that i^K^i produces no work, we obtain: 



dK _ /dW\ 
~dt ^ V~rfr J 



(1) 



(dW\ 



(2) 



+ SE{t) (10) 



where K is the kinetic energy of the polymer and the 
right-hand side of this equation contains the change in 
energy of the U and T sections of the polymer due to the 
mechanical work W'^^\ the viscous drag, T4^(^) = Sh — 
2jK, and the rate of work done by the electric field, Se- 
By writing Eq. (fTU]) we assume that the translocation 
time t{No) for the entire polymer is determined by two 
separate contributions. The first stems from the change 
in free energy of the polymer due to the removal of beads 
from the U side and their addition to the T side. The 
second term arises from shifting the center of mass of 
each blob from the initial position towards (for the U 
part), or away from (for the T part), the pore entrance. 

The simulations reveal that Sh and Se are each inde- 
pendently constant in time to a very good approximation. 
This holds for both the average values over all samples 
(see Fig. [51(a)) and for any polymer realization. The av- 
erage K (not shown) is also approximately constant, as 
the temperature is constant, leading to dK / dt — 0. 

Regarding the mechanical work VF'^^^ simulation data 
show that, except for a short-lived transient at the be- 
ginning and final part of the process, the rate of re- 
moval/addition of beads is linear in time. Similarly, for 
the viscous drag the shift of the center of mass takes 
place at constant velocity, as illustrated in Fig. O by the 

!JJ\ 17^) 

two individual components Sf^ (t) and Sf^ {t) , both of 
which show linear behavior. Since the two contributions 
can be viewed as independent components of the work, 
we separately analyze their effects on the scaling depen- 
dence. 

In order to estimate the work required to shrink/grow 
the two blobs, we now introduce a macroscopic pic- 
ture according to which, for a spherical blob of radius 
Ri, surface Aj and volume Vi the work is dW^^"* = 
PidVj + a-ydAj with the average surface tension and 
Pj — 2a ^/Ri the pressure acting on the blob (Laplace 
equation). Given the anisotropy of the I = U and 



6 




FIG. 4: (Color online) Longitudinal (filled circles) and trans- 
verse (open squares) components of the gyration tensor for (a) 
the U and (b) the T segment with the number of untranslo- 
cated (Ny) and translocated (Nt) beads, respectively. The 
dashed lines show the scaling. 



T segments, the above relation generalizes to dW'-'^'^ = 
J2iJ2f_i=±,\\^i,f^Ri,f^dRi^f_, with A/,^ collecting aU con- 
stants. More explicitly, Ru.A_dRu^± N'^^^'^dNu and 
RT,±dRT.± ^ N^^'^-'^dNr, with i^u - 0.6 and i^t ^ 0.5, 
where the longitudinal components of both U and T seg- 
ments have been neglected in view of their much weaker 
dependence on the number of beads. The rate of work 
on the entire polymer, consisting of the U and T blobs 
with radii given by Eq.®, takes then the form 



dt 



-idr 2 



(1 



with r = NTit)/No. Since dW'^^^dt is constant, in- 
tegration of the above expression, with r £ [0, 1] and 
t G [0, r], leads to the scaling of the total translocation 
time: r ~ Nq'^" . It is worth mentioning that, at variance 
with a previous argument Q, our macroscopic picture 
does not require that each blob be in a state of mechan- 
ical equilibrium. In fact, it is clear from Fig. [3] that 
at the end of translocation the blob is definitely not in 
equilibrium. 

We next consider the viscous drag due to the net mo- 
tion of the blob relative to the fluid, which can be com- 
puted by the global friction experienced by the whole set 

of A^'o monomers Fdmg = J^fj^ili^i ~ ^i)^ where again 
Vi is the bead velocity, Ui is the fluid velocity at the 
monomer location and 7 is the friction coefficient. This 
is best recast in the form of an effective friction coef- 
ficient Fdrag — ~lef fV , whcrc V is the center of mass 
velocity of the blob. As is well known [s^, in the ab- 
sence of hydrodynamic correlations all monomers behave 
independently, so that -feff scales like A^o- On the other 
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FIG. 5: (Color online) (a) Synergy factors with time for the 
hydrodynamic (S^', S^', total Sh) and electric field (Se)- 
(b) Probability distributions of Sh, Se during translocation 
events. Curves are averages over 100 events for TVo — 300. 



hand, when hydrodynamic correlations are included, the 
inner monomers are screened out from the outer ones, so 
that the effective friction is reduced and scales less than 
linearly with length. More precisely, 7e// ~ i? ~ Nq. 
Given that V = dR/dt, the hydrodynamic drag scales 
like iVg^" and A^q", without and with hydrodynamics, 
respectively. Remarkably, the exponent with hydrody- 
namics, 2;/, is exactly the same as the one associated with 
the thermodynamic work W^^\ so that a = 2v in either 
case. In the case without hydrodynamics, however, the 
thermodynamic work and the work due to viscous drag 
exhibit two distinct exponents, 2v ~ 1.2 and 1 -I- z/ ~ 1.6, 
which explains why any attempt to represent the scaling 
through a single exponent is bound to work only on 
a narrow range of values of Nq. Generally, the scaling 

I^Tr(A'o) ^ Nq will be a weighted average of the two, i.e. 

^2z^ < C < 1 + 1^- 

In support of the previous interpretation, we have mea- 
sured the typical values of the rate of change of the blob 
radius Rj, the mean center- mass velocity V and the av- 
erage flow speed U and found that i?/ ~ V/ ~ 5W/. This 
corroborates the idea put forward in this paper according 
to which both mechanisms, the blob shrinking as well as 
the blob shifting processes, must be taken into account 
to provide a complete picture. 



VI. DISCUSSION AND CONCLUSIONS 

Summarizing, we have investigated the process of poly- 
mer translocation through a narrow pore using a multi- 
scale approach which explicitly accounts for the hydro- 
dynamic interactions of the molecule with the surround- 
ing solvent. The translocation time was found to obey a 
power-law dependence on the polymer length, with an ex- 
ponent a — 1.28 ±0.01, in a satisfactory agreement with 
experimental measurements and other computer simula- 
tions. Moreover, our simulations reveal that the coupling 
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of the molecular motion to hydrodynamic correlations 
results in a significant acceleration of the translocation 
process. The scaling behavior observed in the numeri- 
cal simulations has been interpreted by means of a new 
phenomenological model, accounting for the anisotropy 
of both translocated and untranslocated segments. This 
ingredient appears to be crucial to the correct interpre- 
tation of the basic mechanisms behind the physics of the 
polymer translocation, which involves two separate pro- 
cesses, the shrinking and the shifting of the blob. 

Deviations from the mean-field picture occur mainly 
near completion of the process, where the radius of the 
untranslocated segment undergoes an accelerated deple- 
tion: for r{t) > 1/2, the majority of the beads have al- 
ready translocated and entropic forces cooperate with the 
electric field to complete the translocation. Besides vio- 
lating the static scaling at the end of the translocation 
process, entropic forces may lead to more dramatic ef- 
fects, which escape any mean-field description based on 
the translocation coordinate r{t) alone. Such bcyond- 
mean-field-theory effects produce rare retraction events: 



the polymer occasionaly anti-translocatcs after having 
partially passed through the pore. Our simulations reveal 
that retraction events are typically associated with the 
T part entering a low-entropy (hairpin-like) configura- 
tion, which is then subject to a strong entropic pull-back. 
These non-perturbative events depend on the polymer 
length, the initial configuration and the values of other 
parameters (friction constant, temperature and strength 
of the pulling force). They occur at a rate up to 2% 
and do not significantly affect the statistics of the scaling 
exponent. 
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